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ABSTRACT 

We report the discovery of peculiar features in the optical spectrum of 4C+22.25, a flat spectrum 
radio quasar at z = 0.4183 observed in the SDSS and in a dedicated spectroscopic follow-up from 
the Nordic Optical Telescope. The R/3 and Ha lines show broad profiles (FWHM-12,000 km s" 1 ), 
faint fluxes and extreme offsets (Av = 8, 700 ± 1,300 km s _1 ) with respect to the narrow emission 
lines. These features show no significant variation in a time lag of ~ 3.1 yr (rest frame). We rule out 
possible interpretations based on the superposition of two sources or on recoiling black holes, and we 
discuss the virtues and limitations of a massive black hole binary scenario. 
Subject headings: quasars: individual (4C+22.25) 
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1. INTRODUCTION 

4C+22.25 (RA: 10:00:21.8 Dec: +22:33:19 (J2000.0)) 
is a fiat spectrum radio quasa r at z = 0.4183, discove red 
through radio observations bv lMerkeliin et all (11968f). A 
first optical spectrum was collected by ISchmidtl (1974) 
who observed a flat continuum with no significant emis- 
sion line, suggesting that the source is a BL Lac object. 
lHaddad fc Vanderriestl (|l99l re-observed the 4000-6000 
A range and detected a set of bright, marginally resolved 
narrow lines ([Nev] 3 34 6] 342 6 ; [Neni]3 869i 39 68 ; [On] 37 27; 
H7, [Oiii]4363) at z = 0.419. The intensity ratios of 
these lines suggested that the source hosts a S eyfert-like 
Narrow Line (NL) region. iNilsson et al.l (|2003f ) collected 
ground-based high-resolution images of 4C+22.25 as a 
part of a study of blazar host galaxies, and showed that 
the host galaxy is well resolved. Its light profile follows 
a de Vaucouleurs law with scale radius R e = 3.3 ± 0.2" 
(18 ±1 kpc) and apparent magnitude tub, — 18.63 ±0.05. 
Including a fc-correction (0.7 mag assuming a typical El- 
liptical galaxy spectrum at z = 0.419), the inferred lu- 
minosity is Mr = —23.9. 

A companion galaxy is located ~ 6" North- West of 
4C+22.25 (~ 30 kpc at th e redshift of the quasar). 
lHaddad fc Vanderriestl (|1991| ) reported the detection of 
the Call (H) and G-band features in its spectrum, yield- 
ing z — 0.416, and suggested that a gravitational in- 
teraction with the quasar host galaxy may be occurring. 
SDSS photometry reveals also the presence of 12 galaxies 
within a projected separation of 400 kpc and photometric 
redshift consistent with the one of 4C+22.25, suggesting 
that the source may be located in a relatively rich galac- 
tic environment. 
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In this Letter we report the discovery of peculiar broad 
lines in the optical spectrum of 4C+22.25 that is pub- 
licly available fro m the Sloan Digital Sky Survey (SDSS; 
lYork et all 120001) database. Very broad and rather faint 
Mgn, H/3 and Ha lines are observed, all showing a ve- 
locity blueshift of w 8,700 km s _1 with respect to the 
NLs. 

Similar velocity offsets have b een already observed 
in a handful of SDSS quasars: iKomossa. Zhou fc Lul 
(|2008f ) reported a - 2,650 km s" 1 shift between the 
main narrow line system and a second set of narrow 
and broad lines in SDSS J092712.65+294344.0 (here- 
after, J0927), which has been i nterpreted as the sig- 
naturc of a recoiling black hole (|Komossa. Zhou fc Lul 
2008) , of a massive black ho l e binary (BHB; 
i Bogdanovic. Eracleous fc Si gurdsson 120091: iDotti et al 



2009) or the superpos ition of two ob j ects ( Heckma n et al 
2009). Similarly, IShields et all ((20091) found a ~ 
3, 500 km s" 1 shift between narrow and broad lines 
in the spectrum of SPSS J105041 .35+345631. 3 ( here- 
after, J1050). Finally, iBoroson fc Lauerl (pCJOg ) re- 
vealed the presence of a peculiar profile in the broad 
lines of SDSS J15363 6.22+044 127.0, which could b e 
due t o a BHB (Bor oson fc L aucr 2009; Laucr & Boroso 
2009) , a superposit ion of quasars ( Wrobel fc Laorll20CI 



Decarli et al.ll200"9l). or an extreme double-peake d emit 



pe 
ter (jTang fc Grindlavl[20CJ9llChornock et al.ll2010T ). How- 
ever, we show here that most of these interpretations are 
unsuitable for 4C+22.25. 

The structure of the Letter is the following: in Sec- 
tion [5] we analyse the SDSS spectrum and present new 
observations collected at the Nordic Optical Telescope 
(NOT). In Section |3] we discuss possible models to in- 
terpret the peculiar features of this source. Conclusions 
are summarized in Section |4] Throughout the Letter we 
will assume a standard cosmology with Ho = 70 km s _1 
Mpc" 1 , fi m = 0.3 and tt A = 0.7. 

2. THE SPECTROSCOPIC OBSERVATIONS 
2.1. SDSS spectrum 

The SDSS spectrum of 4C+22.25 was collected on Jan- 
uary, 2 nd , 200 6 and was published in th e SDSS Sixth 
Data Release (jAdelman-McCarthv et aLI I2008D . SDSS 
spectra have A/AA ~ 2000 and cover the 3800-9000 A 



Decarli et al. 



range. Uncertainties on wavelength calibration amount 
to 0.05 A, while flux calibration formal errors account to 
5%. The signal-to-noise ratio per pixel at 6400 A is 21. 

Figure Q] shows the SDSS spectrum (top panel) and 
the identification of main emission lines (bottom panel). 
Emission lines were fitted with a double-gaussian pro- 
file, following iDecarli et al.l (|2008l) . Relevant informa- 
tion are provided in Table [1] Typical uncertainties in 
the line FWHM are around 10%. NL peak wavelengths 
have uncertainties of 10 to few hundred km s _1 , depend- 
ing on the line flux (see Table [TJ . The broad line peak 
wavelengths are poorly constrained: for H/3 and Ha, 
we estimate recessional velocity uncertainties of 1,900 
and 1,700 km s^ 1 respectively. The narrow lines are 
marginally resolved. Their mean redshift is (z) = 0.4183. 
The [O in] /[On], [Om]/H/3 and [Nev]/[Nem] flux ra- 
tios confirm the presence of Seyfert-like i onization condi- 
tions in the NL region (see F igure [2] and Hcckmanl ll98Ct 
IHaddad fc Vanderriesdfl99l . 

For the first time, we report the detection of broad 
lines in 4C+22.25. The Ha, H/3 and Mgn lines are 
clearly visible, while broad components of other Balmer 
lines and the iron multiplets are too faint to be detected. 
Both Ha and H/3 are very broad (FWHM-12,000 km 
s _1 , i.e., larger than 9 6% of the qua sars in the huge, 
SDSS-based dataset by Shen et al. 2010) and faint with 
respect to, e.g., the narro w [Qui] lines (see Figure [2j 
only 1.1% of the quasars in Sh en et al~ll2010L has larger 
[O III] /H/3 (broad) values). The most striking property 
of these lines is that they show enourmous blueshifts 
(8,700 ± 1,300 km s" 1 ) with respect to the narrow- 
line system. Similar properties (in terms of fluxes, line 
width and shift) are reported also for the Mgn line, 
but since the peak is close to the range covered by the 
SDSS spectrum, the line characterization is not feasi- 
ble with the available data. We use the line width and 
luminosity of broad H/3 to com pute the mass of the ac - 
tive black hole, following [Vestergaard fc Peterson! (2006): 
Mbh=1 x 10 9 M w . Assuming the bolometric correction 
factor bv lRichards et all (|2006f) for the continuum lumi- 
nosity at 5100 A, this yields L/L^dd = 0.035. 

Few absorption features are also tentatively reported, 
namely the Mgn doublet, the Ca (H) and (K) and the 
NaD line, at a redshift consistent with the narrow emis- 
sion lines. 

2.2. NOT spectrum 

We re-observed 4C+22.25 using the Andalucia Faint 
Object Spectrograph and Camera (ALFOSC) mounted 
on the 2.56m Nordic Optical Telescope (NOT) on June, 
2 nd , 2010, i.e., 1,612 days after the acquisition of 
the SDSS spectrum (1, 137 days in the rest frame of 
the source). Long-slit spectroscopy configuration was 
adopted. Grism # 5 yields a spectral resolution A/AA « 
410 (1.0" slit) in the spectral range 5,500-10,000 A. The 
total integration time (45 min) was split into 3 expo- 
sures to allow an easy cleaning of cosmic rays. Stan- 
dard IRAF tools were used to reduce data. Wavelength 
calibration was performed using Th-Ar arc spectra, and 
cross-checked using the sky emission lines in the sci- 
ence spectra. Wavelength residual rms is 1.3 A. Flux 
calibration was achieved observing a spectrophotometric 
standard star. Galactic extinction was accounted for ac- 
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Fig. 1.— Top panel: SDSS (solid line) and NOT (dotted line) 
spectra of 4C+22.25, shifted to the rest frame assuming z = 0.4183. 
The NOT spectrum is shifted downwards for the sake of clarity. No 
significant difference in the two spectra is reported. Bottom panel: 
the fitted components. Thick lines mark the broad lines, thin lines 
mark the narrow features. The two vertical dotted lines show the 
rest- frame wavelengths of Mgn and Ha. The Mgn line model 
is poorly constrained since the line lays at the edge of the SDSS 
coverage. The velocity shift between the broad and the narrow line 
systems is clearly apparent. 



cording to Sc hlegel Finkbeiner fc Davis! (|1998l ). assuming 
Rv = 3.1. The final spectrum is shown in Figure [I] top 
panel. Its signal-to-noise per pixel at 6400 A is 23. The 
NOT spectrum is in excellent agreement with the SDSS 
observation both in terms of fluxes and peak wavelengths 
of the observed features (see Table [T]). Since the SDSS 
spectrum has a better global quality, we will refer to ve- 
locities and fluxes derived from the SDSS data in the 
following analysis. 

3. DISCUSSION 
3.1. What 4C+22.25 cannot be 

A simple explanation of the two redshift systems ob- 
served in 4C+22.25 would be that the broad and narrow 
emission lines belong to two different objects, superim- 
posed along the line of sight. This scenario is disfavoured 
by the lack of narrow emiss ion lines at the redshift of the 
broad line system (see also lBoroson &: Lauerl [20091 ). We 
estimate that a narrow line as faint as 3.6 x 10 40 erg 
s _1 would be detected at 1-sigma with respect to the 
noise of the SDSS spectrum. This limit corresponds to 
0.017x the flux of the broad component of H/3. From 
Figure [2j it is apparent that the number of quasars with 
[Om]/H/3(broad)< 0.017 is negligible. 

Moreover, in order to get both the sources within 
the fiber aperture of the SDSS, their separation should 
be < 1.5", yielding a solid angle < 5.5 x 10~ 7 deg 2 . 
The number density of Active Galactic Nuclei (AGN) 
at 0.35 < z < 0.45 (i.e., in a velocity space three times 
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TABLE 1 

Summary of the main emission lines as observed in the SDSS and 

NOT spectra of 4C+22.25. (1) Line identification. When both 

broad and narrow components are available, they are marked 

with 'b' and 'n' respectively. (2) Observed peak wavelength. (3) 

Redshift corresponding to the observed peak wavelength. (4) 

Velocity difference with respect to the mean redshift of the 

narrow line system, (z) = 0.4183. Negative values correspond to 

blueshifts. (5) Full Width at Half Maximum of the fitted lines. 

Note that no correction for spectral resolution is applied here. (6) 

Line luminosity. 



Line 


-^obs 


z 


Av 


FWHM 


log -Llinc 




fA] 




[km s _1 ] 
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[erg s" 1 ] 
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(3) 
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(5) 
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SDSS spectrum 










[Nev] 


4745.8 


0.4180 


-60 ± 720 


580 


41.60 


[Nev] 


4860.2 


0.4183 


-10 ±120 


560 


41.87 


[On] 


5288.2 


0.4183 


+10 ±50 


190 


41.96 


[Nein] 


5488.9 


0.4184 


+20 ± 160 


110 


41.53 


H 7 (n) 


6155.4 


0.4178 


-120 ±600 


1100 


41.20 


[Om] 


6192.3 


0.4188 


+110 + 500 


910 


41.36 


Hen 


6647.4 


0.4183 


-10 ±300 


280 


41.08 


H/3(b) 


6679.6 


0.3736 


-9700 ± 1900 


12000 


42.34 


H/3(n) 


6897.4 


0.4184 


+30 ± 130 


570 


41.68 


[Om] 


7035.2 


0.4183 


0±30 


130 


42.25 


[Om] 


7103.0 


0.4183 


-8 ±11 


410 


42.73 


Ha(b) 


9078.8 


0.3830 


-8000 ± 1700 


12700 


43.18 


NOT spectrum 










H/3(b) 


6681.6 


0.3740 


-9400 ± 1600 


13000 


42.32 


[Om] 


7104.2 


0.4185 


+50 ± 30 


770 


42.74 



as large as th e velocity offset obser ved in 4C+22.25) is 
~ 0.37 cleg" 2 ([Schneider et al.ll2010D . Hence, the proba- 
bility of having a random superposition is ~ 2 x 1CP 7 , i.e., 
completely negligible if compared to the number of SDSS 
AGN in this redshift bin (~ 3,300). The probability of 
alignment of two AGN substantially increases if they be- 
long to a common physical structure, e .g., a cluster of 
flaxy. This scenario was proposed by iHeckman et al.l 
)09) to interpret the two redshift systems observed 
in J0927, but subsequent ob servations revealed that no 
signifi cant cluster is present (jDecarli. Reynolds fe Dottil 
2009). The 'superposition in a cluster' argument cannot 
be applied to 4C+22.25, as the velocity difference be- 
tween the two line systems is too high to be attributed 
to the potential well of a single physical structure (see 
the statistical analysis bv lDotti fc Ruszkowski|[2010l ). 

Another scenario suggested to explain the velocity 
shifts between narrow and broad lines observed in J0927 
and J1050 is that the black hole in these quasars is 
recoiling, as a result of the coalescence of a BHB 
(IKomossa. Zhou fc Lul [20081: IShields et all I2009T) . The 
maximum recoil achievable during BH coalescence is < 
4,000 km s" 1 dBakeret al.1 120071: iHerrmann et al.l | 2007 ; 
Campanelli et al. 2007; Schnittman & Buonanno 2007; 
Lousto fc Zlochowerll2009t Ivan Meter et al.l l2010( Pl. As 



a consequence, the recoiling scenario is ruled out for 
4C+22.25. 

8 Note that hydrodynamical and/or purely relativistic ef- 
fects can strongly suppress the kick mag n itude ( Schnittman 
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Fig. 2. — The [O III] /H/3 (narrow) flux ratio plotted as a function 
of the ratio between broad and narrow H/3 components. The hori- 
zontal line marks the separation between Seyfert and LINER ion- 
ization conditions (sec Hcckman 1980). Diagonal lines show the loci 
of various [Offl]/H/3(broad) flux ratios. The big dot marks the po- 
sition of 4C+22.25. For reference, we plot with small gr ey dots the 
flux ratios of the 18,101 quasars from the compilation by Sh en et al.l 
1120100 for which the three lines were detected. 4C+22.25 shows 
strong [O III] and very faint broad H/3 as compared to the average 
quasar population. 



3.2. What 4C+22.25 might be 

A possible alternative is that 4C+22.25 hosts a binary 
black hole. In this picture the primary, more massive 
BH resides at the center of a circumbinary gaseous disk, 
located in the nuclear region of the host galaxy, while 
a secondary black hole orbits around it. Because of its 
motion, the secondary black hole simultaneously accretes 
and prevents the primary one from accreting. The veloc- 
ity shift between narrow and broad lines is then due to 
the Keplerian velocity of the secondary black hole with 
respect to t he barycenter of the binary (for more de- 
tails, see e.g. [B ogdanov ic, Eracleous fc Sigurdssonll2009l ; 
IDotti et al.l I2009f h Assuming circular orbits, the orbital 
period t would be: 



t = 2tt 



GM 2 (sin $ cos <j) 
q{l + q) 2 (Av) s 



(1) 



2004; Boedanovic, Reynolds & Miller 2007; Dotti et al. 2010; 
Volontcri et al. 2010; Kcsdcn et al. 2010). 



where Mi. 2 is the mass of the primary and secondary 
black holes, q = M2/M1, d is the inclination angle of the 
rotational axis with respect to the line of sight, and <j) is 
the orbital phase (defined so that <\> — at the orbital 
node maximizing the blueshift of the broad lines). Simi- 
larly, the separation a between the two black holes would 
be: 

_ GM 2 (sin 1? cos <f>) 2 

a - q(l + q)(Av)Z (Z) 

In order to charachterize the properties of the BHB, we 
therefore need an estimate of Mi and M 2 , which are un- 
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known. Following Dec arli et al.l (|2010D . we use the host 
luminosity to infer the expected mass of Mi, assuming 
iljRH /Mhnst.=0-0015 as obser ved in the Local Universe 
(e.g., iMarconi fc Huntl 120031 ) . For an old host galaxy 
stellar population, we infer Mi = 2 x 10 9 M Q . Assuming 
M-2 = 1 x 10 9 M , as derived in Section [2] we obtain 
separations of 0.04-0.08 pc and orbital periods of 15-35 
years for ■& — 45° — 90° and <f> — 0. On the other hand, 
the velocity shift observed in the SDSS and the NOT 
spectra is unchanged within the uncertainties (~ 2, 000 
km s _1 ). This implies that the period should be > 30 
yr. 

Small but not negligible eccentricities arc expected in 
very massive BHBs, driven by three body interacitons 
with stars. For q ~ 1 (as in the present case), the max- 
imum expected eccentricity is 0.1-0.3, depending on the 
mass of the binary and the steepness of the radial dis - 
tribution of stars in the host galaxy (see Sesana 2010). 
Such small eccentricities do not change significantly our 
estimate^ 

We point out that, at these tiny separations, the broad 
line region is expected to be perturbed. This would 
explain the faintness of the broad lines with respect 
to the narrow lines (se e Figu re [2]). In this case, the 
iVestergaard fc Peterson! (|2006l ) recipe used to estimate 
Mi may not be valid. We therefore adopt a different 
role-of-thumb approach to estimate M%, namely assum- 
ing that the quasar is accreting at 10% of its Eddington 
luminosity. In this case, M 2 « 3x 10 8 M©, a w 0.05-0.1 
pc and t ~ 20 — 60 yr. We conclude that the BHB sce- 
nario is a viable one for 4C+22.25. 

4. CONCLUSIONS 

We present the discovery of extremely peculiar features 
in the optical spectrum of the flat spectrum radio quasar 
4C+22.25. The narrow lines are very bright and reveal 
the presence of a Seyfert-like nucleus. Its broad lines are 
faint and flat (FWHM~12,000 km s" 1 ), and blueshifted 
with respect to the NL of 8, 700 ± 1, 300 km s" 1 . This 
velocity offset between broad and NLs is so high that 
scenarios involving a superposition in a cluster or a re- 



coiling black hole are ruled out at high confidence. The 
probability of a chance superposition of two AGN on cos- 
mological scales is so small that it disfavored, especially 
if coupled with the non-detection of any narrow emis- 
sion line at the redshift of the broad line system. The 
massive black hole binary scenario holds for 4C+22.25, 
but the observation of the target in two different epochs 
separated by 3.1 yr (rest frame) allowed us to set strong 
constraints on the possible orbital configurations. New 
observations with a longer time lag will help clarifying 
if the binary model is correct or not. Moreover, obser- 
vations at higher frequencies, e.g., in the X-rays, would 
help in constraining the mass and Eddington rate of the 
accreting black hole. 

Whether 4C+22.25 is a lone object, or just an ex- 
treme case of a new subclass of AGN, including J0927 
and J1050, is not clear, and demands further investiga- 
tion both from a theoretical and an observational point 
of view. 
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